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Abstract. We consider a finite sub-chain on an interval of the infinite XXX 
model in the ground state. The density matrix for such a subsystem was described 
in our previous works for the model with inhomogeneous spectral parameters. In 
the present paper, we give a compact formula for the physically interesting case 

^sO \ of the homogeneous model. 
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1. Introduction 

The present paper continues the study of correlation functions for integrable spin 
chains launched in ^ E] . In our previous works, we found an exact expression 
without involving integrals for the density matrix of a finite sub-chain of the infinite 
fx. \ XXX, XXZ and XYZ chains in the ground state. More precisely, we treated in- 

homogeneous models in which each site carries an independent spectral parameter. 
The problem of finding a compact form for the answer in the physically important 
homogeneous case remained open. This is exactly the problem which we solve in 
the present paper. We shall consider the simplest case of the XXX model. 
Consider the isotropic Heisenberg antiferromagnet with the Hamiltonian 
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Take a finite sub-chain consisting of sites % = 1, • • • ,n. The density matrix p n for 
this sub-chain in the infinite environment is an operator acting on (C 2 )® n . Its matrix 
elements are given by the ground state average 

Here |vac) is the anti-ferromagnetic ground state, e 3 -,£j = +, — , and (Ef) { signifies 
the matrix unit (5 ae 5bi) a b=+ _ acting on the z-th tensor component. It is important 

to consider a vector h n belonging to ((Q 2 )® 2 ™ instead of the matrix p n acting in 
(C 2 )®" . The vector h n is given by: 

n 
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In the sequel we refer to the tensor components of (c 2 )® 2 ™ by the indices 1, • ■ ■ , n, n, • • 
read from left to right. The main result of our previous papers can be formulated 
as follows: 

(1.1) h n = e n "s n . 

Here s n = FJ Sjj and Sjj signifies the s^-singlet \ (v + <g> u_ — u_ ® v + ) in the tensor 
product of two spaces j,j, being the standard basis of C 2 . The operator Q n 

will be defined later (see (I2.5|) \ It satisfies the condition 

(1.2) fil t]+1 =0. 
So, the series for the exponential (jl.lj) terminates. 

2. Inhomogeneous case 

Let us introduce inhomogeneity parameters Ai, • • • , A n to the corresponding sites 
of the lattice (see pQ for more details). Then the operator fl n (Xi, ■ ■ ■ , X n ) becomes 
dependent on these parameters. We define this operator following our previous 
works, but we shall slightly change the notation. 

Let {S a }\ =l be a basis of sl 2 satisfying [S a , Sb] = 2ie a b c S c . Following jl] define the 
L-operator which belongs to t/(s[ 2 ) <S> End (C 2 ): 

(2.1) L(A) = 4^ (0) (A), L^(\) = \ + 1 - + 1 -j2S a ®- a , 

A + 2 1 1 a=l 

where d is related to the Casimir operator as Ylt=i = d 2 — 1, and 

P& d ) = - -71 H ll ) \ l ( ■ 

r(l-f-|)r(i-f + |) 

In this normalization we have the unitarity and crossing symmetry in the form 

L(X)L(-X) = 1, a 2 (L(A))* a 2 = -L(-X - 1). 

We shall consider tensor products of several spaces C 2 . In that case the index i in 
Lj(A) denotes the tensor component as usual. In what follows the function p(X,d) 
always comes in the combination 

p(X,d) p(X-l,d) 1 

A + f-1 ""A 2 -f 

so the T-functions will never really appear. 

We shall also use the ordinary 4x4 i?-matrix obtained as the image of the L- 
operator (J2.1)) in the 2-dimensional representation of [/(s^). When acting in tensor 
product of two spaces it will be denoted by Rij(X). We denote the corresponding 
factor by 

p(X) := p(X, 2) . 

Let us explain the results of the papers [TJ |2] in the setting of the XXX model. 
First, introduce the operator 

TW(A; A 2 , . . . , A n ) := L- 2 {X - A 2 - 1) • • • L n (X - X n - l)L n (X - X n ) ■ ■ ■ L 2 (X - A 2 ) . 
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Notice that in this product the sites 1, 1 are omitted. 

In the paper [Tj, we discussed in detail the linear functional on [/(s^) called Tr\. 
Denote by Wd the irreducible d- dimensional representation of Uisi-z). By definition, 
the map Tr A : £/(s[ 2 ) — > C[A] associates with each A 6 U(sl 2 ) a unique polynomial 
Tr x (A) of A, such that Tr d(A) = tr C dZUd(A) holds for any positive integer d. Some 
of its main properties are 

Tr x (AB)=Tr x (BA), 
Tr x (C£S 2 a )A) = (\ 2 -l)Tr x (A), 



Tv(e ts3 ) 



a=l 

sinh tX 



sinht 

In the present paper we shall also use the linear functional 

Tr A ....... v , : U(sl 2 ) 9k ^C[\ lr .. ,Xk) 

defined by 

k 

Tr^...,^ <g> • • • ® A fc ) = JjTr^.^). 

j'=i 

The main ingredient of our construction is the operator 
(2.2)X n (A l5 --- .An)!,...,^...,! := (-l) n - 1 res A=A2 Tr Al _ A (rW(^^))p 1) 2T- I T- 5 , 

where Pjj is the permutation and O 5 "^ = (1 — Pjj)/2 is the skew-symmetriser. Define 
further 

ClF(A 1; • • • , A n ) = -4w(A iJ )xW)(A 1 , • • • , A n ) , 
X n ^(A 1 ,---,A„) = RP(A 1 ,---,A n ) 

X X n (Aj, Aj, Ai, • • • , Ai, • • • , Xj, ■ ■ ■ , An) i j )1 ...7...^... n ^ ) ...|...f... I j-lRn J (Al, • • • , A n ), 



where 



and 



R^(Ai,..-,An) 

:= -Ri,i-l(Ai,i-i) • • • Ri,i(Xi t i) ■ R—^Xi-ij) ■ ■ ■ i?ij(Ai 5 j) 

x-Rj^t^A^ij) ■ • ■ Pi+rjiAj+ij) • #i=ij(Ai-i,j) • ■ ■ -Rij(Aij) 
K^(A 1; ---,A n ) 

x -Rj,n(Aj )n ) • • • i2jj+i(Ajj + i) • -Rj J _i(Aj ) 

u(X) = — logp(A)+ 1 



dA torv 7 2(A 2 - 1) ' 
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The definition of differs from the one used in [21 E] 1 by the second product of 
i?-matrices, but the final formula (jl.lj) remains unaltered by this modification. 

The operators fin possess a number of properties the most important among 
which are 

(2.3) [n^>(Ai,---,A n ),fi^(Ai,---,A n )] =0, 

(2.4) n£0(A a , • • • , A^n^CAi, • • • , A n ) = if {i,j} n{k,l}^0. 

As a function of Ai, • • • , A n , f2„ (Aj, • • • , A„) is meromorphic. All poles are simple 
and located at Aj, Xj = A; and X i} \j = Xi ± 1 for I ^ i,j (which are due respectively 
to X^(Ai, ■ • ■ , A n )), and to the i?-matrices) and Ajj G Z\{0} (due to u(X{j)). We 
have 

h n (Xx,--- ,An) = e n ^ Al '-' A »)s n , 

where 

(2.5) n n (X ir -- ,\ n )= J2 ^ ,i} (Ai,--- ,A n ). 

l<i<j<n 

The properties ([2.3)1 . (|2.4|) guarantee the nilpotency (|1.2j) . 



3. Homogeneous case 

Our goal is to obtain the homogeneous limit Ai = • ■ ■ = A n = 0. In the original 
formula (|2.5|) . this problem is very complicated: the singularities on the diagonal 
Aj = Aj are present in every term of ()2.5|) . Although these poles are absent in the 
sum itself, it is technically difficult to explicitly carry through the cancellation and 
obtain the final answer. 

So, we need to rewrite the formula for Q n (X\, • • • , A n ) in such a way that taking 
the homogeneous limit is easier. To this end, let us write first of all another formula 
for X n (Ai,-- - ,A n ). 

Denote by zu x the A-dimensional irreducible representation. We use only the fact 
that the Casimir element reduces to A 2 — 1, and hence the following computation 
makes sense for non-integer A as well [1 . 

Notice that 

^ + (A) = i«7 A (L<°>(A/2)), P-(A) = -~^a(4 0) (-A/2)) 



See eq.(12.1) in We have also used the fact that, in the notation there, Jva(Ai,2/2)£2(^i,2/2— 
1)T~2 = inside the trace and P2,2 s i,2 s I,2 ° n-2^n = Pi^ii^i' 
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are orthogonal projectors. Consider now X n (X 1: • ■ • , A n ) . Using the formula Pi^i {$2 2 
t y^ 2 t y^^Pi 2 1 the definition of the L-operator and the crossing-symmetry, one finds 

-m, 2 (4 0) (^)4 0) (^)) = - ^ (4 0) (^)4 0) (^ - d) 

= - (a 1>2 - 1)^, 2 (4 0) (^)) ^2, 

(4 0) (^ - d4 0) (^ - d) ^ = - (4 0) (^ - D4 0) (^)) ^ 

= (Al,2 + lW. a (^(^f 1 )) ^2 

= "(Ay + IK, ^2- 

Now it is easy to see that 
(3.1) 
X n (X±, • • • , A n ) 

= (-l) n - 1 res /11=Al res M2=A2 - Tr Ml2 (T„(^-^; Ai, • • • , A^P^O^O^ 

where T n (A) is the complete monodromy matrix: 

T n (A; A l5 . . . , A n ) = L\{\ — Ai — 1) • • • L S (A - A n - l)L n (A - A n ) ■■■L 1 (X-X 1 ). 
Using the Yang-Baxter equation one finds the following formula for 

X$*X\ U • • • , A n ) = (-l^res^n*^ -^l^Tr^ (r n {^±^- X u ■ ■ ■ , A n )) 

x R^ } (Ai, • • • , A„)P iJ «- J Ri i ^(A 1 , • • • , A n ). 
By a straightforward computation one finds 

Ri^(Ai, • • • , A, ; )/>.;j' ( ;j' ;; I'i;-sA l . • • • , A n ) 

= 2^ r2 ' 2 C^»»(A»' ' • • ) \i) <H> T n (Xj] Ai, • • • , A n ) • 5 ) . 

Here the skew-symmetriser 7~ acts on the auxiliary space C 2 <E> C 2 . Notice that 
Tr2,2 (T n (iii] Ai, • • • , A n ) <S> T n (/i2] Xi, • • • , A n ) • CP~) is actually symmetric with respect 
to /ii, fi2 due to the relation 

Obviously, the formula for Q n (Xi, • • • , A n ) can be rewritten now as 
(3.2) n n (K ■ ■ ■ , A„) = t^lj j gig 2 . W(M Tr m , 2 (T n (^; Ai, • • • , A„)) 

x Tr 2j2 (T„Gui; Ai, • • • , A n ) <g> T n (/i 2 ; Ai, • • • , A n ) • P(/ii, 2 )) , 
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where 

^1,2 1 

the contours of integration encircle the poles [i\ = Xj, fi2 — Xj for j — 1, ■ ■ ■ , n. 
The great advantage of this formula is that it allows to take the homogeneous limit 
Xj = 0. In the next formula we write Q n for Q n (0, ■ ■ ■ , 0), etc.. 

(3-3) n n = —JJ—— u{jMlfl )Tt^T n (— r -)j 

x Tr 2i2 (r n (^i) ® T ft (/i 2 ) • , 

where the integrals are taken around /i« = 0. 

Formulas (|3.2j) . (|3.3|) are the main results of the present paper. Let us discuss 
them briefly. 

First we note that the integrand of (|3.3j) has a pole of order n at fij — 0. By 
evaluating residues, Q n becomes a linear combination of the Taylor coefficients of 
a; (A) with j < n, given explicitly by 

^ (A) " 2 (A 2 - 1) = 2 ( bg2 + J^ U2k + 1)A2fc ) ' 

Here Ca( s ) = (1 — 2 1_s )C(s), C( s ) denoting the Riemann zeta function. This settles 
the conjecture of jSHH] which states that any correlation function of the XXX model 
can be written as a polynomial of log 2 and C(3), C(5), • • ■ with rational coefficients. 

Second, formula (|3.3|) may open up a way for studying the large-distance limit. 
Also it would be very interesting to see if it helps for the investigation of the limit to 
continuous field theory. We hope to return to these problems as well as the extension 
to the XXZ and XYZ cases in future publications. 
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